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Abstract 

A narrow peak in the tt^tt" mass distribution was seen by the Rome- 
Syracuse Collaboration in pn — >■ 27r+37r~ annihilation at rest in 1970. It 
was ignored for 40 years. The reanalysis of this peak finds that it has 
X^ , the mass 757.4 ± 2.8stat ± 1.2sys MeV/c'^ and a width consistent with 

^s^ ■ the experimental resolution. The evidence of the peak is 5.2 standard 

deviations. The peak is generated in (1.03 ± 0.21stat ± 0.21sys)% of the 
pn annihilations at rest. No spin analysis is possible with the statistics of 
the experiment but there are arguments suggesting that it has J^ — 0^. 



1 Introduction 



In 1970, the Rome-Syracuse Collaboration (RSC) studied the branching ratio 
C^ ' of the decay w — >■ tt+tt" [T] using the data of the pn annihilations at rest 

collected in the 30" BNL bubble chamber. The analysis measured the upper 
limit 4.3% at 95% confidence level and found an unexpected result: the tt+tt" 
mass distribution of 1496 annihilations at rest 

pn -^ 27r+37r" (1) 

had a narrow peak at about 755 MeV/c^. This distribution is shown in Fig. la. 
A x^ fit of this distribution found that the peak had a significance of about 
4.5 standard deviations (SD) and a width lower than the experimental resolu- 
tion. No relation was found between the tt'^tt" and other angular and mass 
distributions. These facts suggested to the RSC that the peak was generated 
by a fluctuation. 



c 250 - 



■-D 200 - 



150 



100 



50 




250 500 750 1000 1250 
tt^tt" moss (MeV/c^_ 



400 600 800 1000 
TT^Ti" moss (MeV/c^) 



Figure 1: (a) The tt+tt" mass distibution of the pn — > 27r+37r^ annihilations 
at rest measured by the Rome-Syracuse Collaboration. The dash-dotted line 
is the phase space prediction; the dashed line is the prediction of the analysis 
of Ref. p]; the solid line is the prediction of the fit F reported in Table I; the 
dotted line is the polinomial background in the same fit. (b) The tt+tt" mass 
distibution of the np — )■ 37r+27r^ annihilations at low momenta measured by the 
OBELIX Collaboration; the crosses are the Rome-Syracuse data renormalized 
to the OBELIX combinations; the solid line is the convolution of the Rome- 
Syracuse data with a Gaussian having the resolution FWHM = 37.0 MeV/c^. 



The properties of the annihilations ([T]) were undestood at the beginning of 
1990s [21 13]. A reanalysis of the RSC data proved that this reaction is dom- 
inated by the channel pn — >■ /o(1370)7r~ followed by the /o(1370) decay into 
p(770)°p(770)*' and S^S^, where S^ indicates the tt+tt" / = S-wave interac- 
tion. The prediction of Ref. [3] are shown by the dashed curves in Fig. la. It 
fits very well the experimental data, with the exception of the interval 700-800 
MeV/c^. The parametrization of Sw used for obtaining this curve is shown in 
Fig. 2. 

A year after the communication of the preliminary results of this analysis 
[2] , the OBELIX Collaboration presented the preliminary results of the analysis 
of the charge conjugate annihilation at low momenta [4j 



np —^ 37r^27r 



(2) 



The 7r+7r~ mass distribution of OBELIX is shown in Fig. lb, togheter with 
the RSC distribution normalized to the OBELIX combinations (crosses). This 
figure shows that the tt+tt" distributions of both experiments are in agreement 
with the only exception of the peak at 755 MeV/c^ that is not shown by the 
OBELIX data. At that time, the absence of the peak in the OBELIX data 
confirmed the opinion that the RSC peak was a fluctuation. 

Recently, Troyan et al. claimed to have observed several narrow tt+tt^ peaks 



in the reaction np 



npn^i: 



with neutrons of 5.2 GeV/c fS^. One of these 
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Figure 2: The parametrization RPOH4/5 of the tt+tt" / = S-wave interaction 
used in Ref. [3]. The dotted hue shows the K^K^ threshold, (a) Phase-shifts, 
(b) Inelasticities. 

peak has a mass close to 755 MeV/c^. The coincidence of the mass of the RSC 
peak with that of Troyan et al. suggested to reanalyse the RSC data. 



2 Reanalysis of the Rome-Syracuse distribution 

The 7r~''7r~ mas distribution of the RSC has been fitted in the mass interval 
600 - 900 MeV/c^ with the parametrization 



D{m) ^P:i{m)+NS{m), 



(3) 



where Pz{m) is a third degree polinomial, S{m) is the signal parametrization 
normalized to one, and N is the number of combinations in the peak. 

The fits have been carried out either with the x^ or with the unbinned 
maximum likelihoocl^ methods. Three signal parametrizations have been used: 

(a) Gaussian with the resolution a free. 

(b) Gaussian with the resolution fixed at cj = 12.74 MeV/c^. (It correspond 
to the FWHM = 30 MeV/c^ estimated by the RSC in Ref. [T].) 



(c) The double Gaussian 



S{m) 



2^? 



27rtTi 



/27rf72 



(Af-m)^ 
2^ 



(4) 



^ The maximum likelihood fits were carried out by using the known values of the tt+tt" 
masses, which were maintained under the support of punched cards and afterwards copied 
manually. 



Table 1: Results of the x^ fits of the distribution shown in Fig. la in the interval 
600 < m < 900 MeV/c^ Fit A is the fit with a free. Fit B is the fit with a 
fixed at 12.74 MeV/c^. Fit C is the fit with the resolution function ^. Xo is 
the x^ of the fit with a third degree polynomial without the signal. The last 
row reports the evaluation of the peak significance. 



Fit 


A 


B 


C 


N (combinations) 


152 ± 33 


188 ±37 


220 ±45 


M (MeV/c2) 


755.4 ±1.9 


756.3 ±2.5 


756.3 ±2.7 


a (MeV/c^) 


7.7±1.7 


12.74 




X^ 


14.7 


19.6 


20.5 


d.o.f 


24 


25 


25 


xVd.o.f. 


0.61 


0.78 


0.82 


X^ 


44.2 


44.2 


44.2 


7V/A7V 


4.6 


5.1 


4.9 



Table 2: Results of the unbinned maximum likelihood fits in the interval 600 < 
m < 900 MeV/c^ Fit D is the fit with a free. Fit E is the fit with a = 12.74 
MeV/c^. Fit F is the fit with the resolution function (j4]). The last rows report 
the evaluation of the peak significance. 



Fit 


D 


E 


F 


N (combinations) 


147 ±35 


186 ±36 


224 ± 43 


M (MeV/c2) 


756. 7 ±1.9 


757.4 ±2.6 


757.4 ±2.8 


a (MeV/c2) 


8.1 ±2.1 


12.74 




N/AN 


4.2 


5.2 


5.2 



with /i = 0.7483, /2 = 0.2517, ai = 13.24 MeV/c^, and (72 = 32.45 
MeV/c^. These number have been obtained by fitting the error distribu- 
tion of the Rome subsample. (665 events.) This function has FWHM = 
33.6 MeV/c2. 

The results of the fits with the x^ ^^nd with the maximum likelihood methods 
are reported respectively in Table I and in Table II. These fits prove that the 
significance of the peak is about 5 SD when the resolution is fixed using the 
signals (b) and (c). In particular, the fits with the maximum likelihood and the 
parametrizations (b) and (c) have both the significance of 5.2 SD. They prove 
also that the evidence of the peak was not understood in 1970s because nobody 
of the RSC realised that the use of the known experimental resolution would 
have improved the significance above five SD! 

The best fit is F, because it takes into account the tails of the error distri- 
bution and is made with the maximum likelihood method. Figure la shows the 
prediction of this fit (solid lines) and of the polynomial background of the same 



fit (dotted line). 

3 Discussion of the difference between the Rome- 
Syracuse and the OBELIX data 

The absence of the peak in the OBELIX tt+tt^ mass distributions can have 
two explanations: (i) the RSC peak was a fluctuation; (ii) the OBELIX Collab- 
oration did not see the peak because the detector had a resolution poorer than 
our bubble chamber experiment. 

The explanation (i) seems improbable because the unbinned maximum like- 
lihood fits with the signals parametrizations (b) and (c) have the significance 
5.2 SD. The explanation (ii) is possible because the OBELIX detector measured 
the charged pion tracks after they had traveled through the liquid hydrogen of 
the target and the material of the cryogenic vessel. 

An indication that the OBELIX Collaboration could not have seen the peak 
because a poorer resolution is given by smoothing the RSC distribution. The 
solid line in Fig. lb shows the convolution obtained by substituting a Gaussian 
having FWHM = 37.0 MeV/c^ to all the tt+tt" masses of the RSC sample. It 
estimates the tt^tt" mass distribution if the RSC data had had the resolution 
FWHM = 50 MeV/c^ instead of 33.6 MeV/c^. (In fact, 50^ - 33.6^ = 37.0^.) 
This curve does not show any peak and is perfectly compatible with the OBELIX 
distribution. 

There is a more convincing argument that supports the explanation (ii). The 
TT+TT ~ mass distribution is given by 

D{m) = I \A{m;{x})\^d<^{{x}) = R{m)^{m) 

where A{m; {x}) is the amplitude of reaction ([T]) that depends on the tt+tt" 
mass m and on a set {x} of other kinematical variables, d^{{x}) is the element 
of the phase space volume, ^(m) is the phase space distribution, and R{m) 
is the value of |A(to;{x})P mediated on the variables {x}. R{m) is the ratio 
between the experimental distribution and the phase space (REDPS) 

J\A{m;{x})\^dmx}) _D{m) 
^""^ Jd^lx}) $(to)' 

Figure 3a shows this ratio for the RSC data. It has two maxima: one at 
about 755 MeV/c^, and the other at about 1250 MeV/c^. The first maximum is 
generated by the peak under study, the second is probably due to the production 
of the /2(1270) meson. This ratio shows also an ankle at the KK threshold and 
does not show any evidence for the debated cr(600) meson. Fig 3b shows the 
REDPS of the OBELIX data. It has the same behaviour shown by the REDPS 
of the RSC, and confirm the peaking at about 755 MeV/c^. But its maximum 
is wider than that of the RSC and is not reproduced by the prediction of fit F. 
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Figure 3: The ratios between the tt+tt" mass distibution and the phase space 
predictions, (a) Ratio of the Rome-Syracuse data, (b) Ratio of the OBELIX 
data. In both diagrams the dashed hue is the prediction of the fit made in Ref. 
^ , the sohd hue is the prediction of the unbinned maximum hkehhood fit of the 
RSC data with the formulae ^ and (j4]), and the dotted fine is the pohnomial 
background estimated by the same fit. 



4 Quantum numbers of the peak 

The REDPS distributions prove that there is an amphtude peaking at about 
755 MeV/c^. This peaking cannot be generated by a reflection because no 
significative correlations were found between the tt+tt" peak and other mass 
and angular distributions. The low mass suggests that it cannot have J > 1. 
Furthermore, the REDPS peak is too narrow for being produced by the p{770)'^ 
resonance and its mass is 25 MeV/c^ far from that of the w meson. Lastly, the 
peak cannot be due to the p(770)" — uj interference. In fact, the annihilations 
(HI are dominated by the S-wave initial states (SJ 1^. Because the G-parity 
conservation, the channel pn — >■ p(770)°7r+27r~ is generated by the ^^o state, 
while the channel pn -^ a;7r+27r~ is generated by the "^Si state. These states 
cannot interfere because they are not coherent. 

Therefore, the best hypothesis is that the peak is generated by the J^ = 
O"*" interaction. It could be generated by an unknown scalar meson or could 
be a property of the Sw interaction that was not reproduced in the previous 
analysis 2, 3 . In fact, the phase-shits parametrization shown in Fig. 2a is 
below 90° tiU up 862 MeV/c^. This hypothesis sugests that the peak could 
be reproduced by forcing the phase shift parametrization to pass at 90° at 
757.4 MeV/^. Unfortunately, the events of the RSC are too low for allowing a 
convincing study of the effects of the variations of the / = S-wave amplitude. 



5 Conclusions 

The fit F measures the foUowing parameters for the narrow tt+tt^ peak shown 
in Fig. la 

M = 757.4 ±2.8,tat±1.2sysMeV/c2, 
r < 30 MeV/c2, 
N = 224 ± 43st,at ± 45sys combinations, (5) 

where the systematic error has been evaluated by using the mean of the differ- 
ence between the parameters found in fit F and those measured in the other five 
fits. The mass M is close to the values measured by Troyan et al. [5]. 

There is no evidence for double peak production in the same event. There- 
fore, one evaluates from ^ that the fraction of the channel P(757)7r"'"27r~ in 
the pn -^ 27r*37r^ is 

N 

— - = (15.0±2.9stat±3.0,ys)%. 
1496 

Since the frequency of the annihilation ([T|) at rest is F{pn — > 27r+37r^) = (6.9 ± 
0.5)% [3], the frequency of the peak is 



sysy 



70. 



F[P(757)7r+27r-] - ^ F{pn -^ 27r+37r-) - (1.03 ± 0.2l3tat ± 0.21 
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